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ABSTRACT 

in the paper, the theory and experimental measurements on the aerodynamic decay of a wake from high per- 
fcrmance vertical axis wind turbine will be discussed. In the initial experimental study, the wake 
downstream of a model Darrieus rotor, 28cm diameter and a height of 45.5cm, was measured in the Univer- 
sity Bouridary Layer Wind Tunnel. The wind turbine was run at the design tip speed ratio of 5.5. It was 
found that the wake decayed at a slower rate with distance downstream of the turbine, than a wake from 
a screen with similar troposkein shape and drag force characteristics as the Darrieus rotor. The ini- 
tial wind tunnel results indicated that the vertical axis wind turbines should be spaced at least forty 
diameters apart to avoid mutual power depreciation greater than ten per cent. 


INTRODUCTION 

An .lmi>ortant aspect concerning the wake of a 
Darrieus rotor, is the interaction of such tur- 
bines in widespread arrays. For the production 
of cheap electrical power from the wind, it has 
been proposed that large wind turbines could be 
placed in ‘clusters’ or ’arrays’. An initial 
study of this problem by R.T. Templin C1974) , 
concluded that if the spacing between individual 
power units was more than about 30 rotor diame- 
ters, the power reduction per machine would not 
exceed 5 to 10 per cent. However, if the spacing 
is further reduced, there is a rapid increase in 
the wake interference. The report treated the 
rotors as an array of large-scale roughness ele- 
ments added to the terrain surface that was al- 
ready rough, and hence, the individual wakes of 
the turbine rotors were 'smeared out* by the 

turbulence she ar layer, well within 
the average spacing between adjacent machines. 
Since this first paper bh the interaction of 
windmills in widespread arrays, several other 
studies have been made, notably by Craford (1975), 
Llssaman (1977) , Newman (1977), Built jes (1978), 
Bragg and Schmidt (1978) , Faxen (1978) , and 
Built jes and Smit (1978) . 

Tn general, the development of the vortex wake 
downstream of a Darrieus rotor in the ambient 
bound ary layer, may be divided into three re- 
gions. The first region represents the formation 
and rolling up of the concentrated vortex sheet 
shed from the rotating aerofoil blades, and then 
the decay of the coherent wake by vorticity am- 
plification and molecular and turbulent diffusion 
and, finally, a turbulent wake uncorrelated with 
tlie initial disturbance. For a two-bladed 
Darrieus rotor when the vertical plane of the 
rotor is normal to the flow direction, the lift 
force on the blades is zero. For the time period 
for the rotor to rotate half a revolution, a 
closed vortex system is shed from the blades and 
is convected downstream, changing shape due to 
the mutual interaction, the effect of the am- 
turbulent boundary layer, and the ground 
effect. In the case, of the two-bladed Darrieus 
rotor, a 'packet' of vorticity is shed periodi- 
cally . 

Hie dominant factors involved in the description 
of the motion within the Darrieus rotor wake, 
are the initial properties of the wake, which 
governs the rate at which the wake entrains the 


surrounding atmosphere, and the rate and manner 
in which the entrained fluid is mixed into the 
wake. For example, it has been well illustrated 
(see Olsen, Goldburg and Rogers) that the two 
trailing vortices in the wake of an aircraft wing 
eventually break down to closed vortex configura- 
tions, The wake of the Darrieus rotor starts 
with a closed vortex structure. On the other 
hand, the wake behind a wire screen is already 
very turbulent, and shows very little evidence 
of any initial coherent vortex structure. It can 
be hypothesized therefore, that the wake of a 
horizontal axis wind turbine, which can be com- 
pared to that of an aircraft wing with the vorti- 
ces forming a helix pattern, decays slower with 
distance downstream than the wake from a vertical 
axis wind turbine, when both turbines are tested 
at the same power coefficient and in the same 
turbulent boundary layer. This fact is important 
when discussing suitable spacing of turbines. 

In the report, a vortex model of the near wake 
of the Darrieus rotor, will first be discussed, 
and finally, results of a study made on the simu- 
lation of the Darrieus wake by arrays of Tropo- 
skein shaped grids will be presented, 

NEAR WAKE OF A DARRIEUS TURBINE - VORTEX MODEL 

In the analysis, the Darrieus rotor troposkein 
shape was approximated by straight line elements. 
Since most of the torque is generated by the . 
’equatorial’ element, the relevant theory will 
be discussed considering this element. In Fig, 

1, a plan view of a typical straight bladed ele- 
ment is shown. The ambient airflow approaches 
from the right, and at the blade, the total velo- 
city is the vector sum of the blade rotational 
speed and the free stream velocity corrected by 
a single inflow or blockage factor. The inflow 
factor takes into account the linear momentum 
change due to extraction of energy by the tur- 
bine. 

Previous theory by Wilson and Lissaman (1974) , 
and also by Base and Russell (1979) showed that 
the inflow factor (a) varied as the turbine soli- 
dity (cf) and Tip Speed Ratio (Tsr) r so that it 
could be suitably represented by the following 
expression: 

a = a (Tsr) I Sin 9| (For 4<Tsr<6) (1) 

Where 6 is the actual blade angle . By consider- 
ing the geometry in Figure 1, then: 
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Tan a == 


U Sin e 


( 2 ) 


(U Cos 8 + Ut) 

Defining the Tip Speed Ratio and inflow factor by: 

m 


and U — Uoo(l”a) 


llien : Tan a 


(l-a)Sin 8 
(l-a)Cos 9 + X TsR 


(3) 


For the equatorial element X = 1, and also siib- 
stituting Equ. (1) into Equ, (3) , then the follow- 
ing equation is obtained for the angle of inci- 
dence of the aerofoil: 



The aerofoil will experience a lift force whose 
magnitude will depend on the incidence of the to- 
tal velocity vector and this lift force will be 
related to the circulation (F) developed about 
the aerofoil. 


According to Lanchester*s wing theory, line vor- 
tices may be considered to exist within the aero- 
foil, These are the so-called 'bound* vortices. 
However, the bound vortices for a finite wing 
cannot end in space, but become free trailing 
vortices ending downstream by a transverse vortex 
or the starting vortex. In plan view, the vortex 
array, in its simplest shape, is a rectangle with 
the trailing vortices linking the bound vortex to 
the starting vortex. 


The concept of a line vortex is derived from a 
vortex tube by making the area of cross section 
of a vortex tube very small, while the circula- 
tion strength *F* remains unaltered, Glauert 
(1946) discussed the line vortex and also the 
equation for the induced velocity at a point P 
of an element of a line vortex at a point Q, which 
is given by: 


dq 


Fds 

TUP 


Sin 0 


(5) 


where F is the circulation defined by J— * — the 
line integral of the velocity around a closed 
curve *c* and for a potential vortex is constant. 
This is the Biot-Savaart Law. The diagram below 
shows an element of the line vortex where 'r* is 
the distance between the field point P and the 
element at Q, and 8 is the angle between the 
direction of the element and the line joining the 
element to the point P. The velocity dq is normal 
to the plane containing r and ds, and its 'sense' 
is the same as that of the circulation *F* of the 
elemental vortex. 



In general, the vorticity has the same geometri- 
cal properties as the velocity field. The 
streamlines of an incompressible fluid therefore, 
correspond to vortex lines of the rotational flow 
of which the direction at every point is that of 
the vorticity vector. An element 'ds' of a line 
vortex therefore, cannot exist independently and 
the Biot-Savaart Law, for an element of vortex, 
should only be used as "building bricks" for 
more complicated closed vortex loops. 

Another expression for the "elemental vortex" is 
given by : 

aq = ^ (6) 

where 'h* is the length of the normal from the 
field point (P) to a line passing through the 
vortex line as shown in the Figure. 

In the vortex model of a straight bladed vertical 
axis wind turbine, the aerofoil blade was as- 
sumed to start at 6 = 0 as shown in Figure 1, 
rotating at a speed of *^2' in a free stream of 
velocity 'Uoo' to achieve a particular Tip Speed 
Ratio. The lift coefficient of the aerofoil at 
this angle (6=0) is zero, and hence, the first 
shed vortex from the aerofoil has zero circula- 
tion. As the aerofoil rotates, by considering 
the angle of the local flow incidence to the 
aerofoil derived in Section 3 and the lift-incl- 
dence curve for the aerofoil, the circulation 
of the shed vortex may be determined as follows. 
The circulation is related to the aerodynamic 
lift on the aerofoil element by the Magnus 
effect equation and; 

Lift Force = pUooF dz 
^ 0 

Hence the lift force is given by the expression 
(assuming no tip losses) 

L = pUc»Fh 

The lift coefficient is defined by the expression 


Cl = 


L 

•|pUx^ ch 


Where Ut is the total velocity of flow approach- 
ing the aerofoil and is the vector sum of the 
free stream velocity and the aerofoil blade ro- 
tation with allowance for the blockage factor. 

It can be shown that the aerofoil circulation is 
related to the aerofoil lift coefficient by the 
following expression: 

r = CLcUoo f uT] 

2 (Ua>J 

The variation of lift coefficient with incidence 
(a) is also known for the particular aerofoil 
section. Finally by previous theory (Equ. 4) , 
the angle of incidence is related to the blade 
angle (9) , solidity (a) , and Tip Speed Ratio 
(Tsr) . 

and also: U-f = (1-0 (Tsr) | Sin9 |)Uoo SinB (7) 

sina 


Figure 2 shows the initial rectangular vortex 
pattern and Figure 3 shows the complete vortex 
model. Further studies on this vortex model are 
currently in progress. 


52 



TESTS TO DETERMINE THE FAR FIELD WAKE 

Tr^e main objective of this experiment was to 
model the interaction that occurs within a wide-* 
3pfe^~t^6-dimensional array of Darrieus wind 
turbines. Figure 4 shows a sketch of a two- 
~5Taded turbine and some of the essential dimen- 
sions, The two-bladed turbine blades consisted 
of extruded aluroinum airfoil section NACA 0012 
wnich had the necessary high-lift, low drag char- 
acteristics, and were curved in the shape of a 
troposkein cantenary. The centre column served 
a| a punting point for the blades and also 
housed two low-friction roller bearings. This “ 
whole assei^ly was then mounted on a centre 
shaft whic h was braced by wires at the top to 
reduce “turbine vibrations. Two vertical axis 
-tu^xTies“used ig the study were geometrically 
“identical#, being 45, 5cm high, 28.0cm diameter, 
and with an aero^foil chord of 2.5cm. This gave a 
sclldity which is defined as the ratio of blade 
area to wind turbine area of ,36. Since the 
THOdels used were operating with little load, the 
Tip Speed Ratio measured was approximately 5,5, 
The power coefficient for this Tip Speed Ratio 
and solidity was approximately 0.4. An important 
problem encountered when working with small 
TRcdels, is that bearing friction will often be 
greater than the power produced, and therefore, 
tlie windEdll will not run. The Darrieus turbines 
used in this study were approaching the lower 
limit of geometric size. It is most difficult 
tc constiruct working models of 2,5cm diameter, 
which would be required in the turbine array 
studies . 

This i>articular experiment was an attempt to 
determine the special requirements for adjacent 
turbines, while allowing for a 10% drop in power. 
This 10% power loss applies to the last windmill 
of the array, so that if there were 10 rows of 
windmills, the windmill in the last row would be 
producing 10% less power than the windmill in 
the first row. In order to determine the power 
loss in a widespread array, it was necessary to 
Hiodel a large array, and then obtain velocity 
measuirements downstream from the array. The 
power extracted by the turbine, is proportional 
to the velocity cubed, and so by obtaining the 
drop in velocity caused by the wind turbines, 
it was possible to calculate the percentage 
power loss. 

In order to perfoim the experiment in the Uni- 
versity Boundary Layer Wind Tunnel, with dimen- 
sions 2.4m X 1.8m x 24.0m, it was necessary to 
use vtjry small models in order to obtain an 
array of acceptable size. It was not possible 
to construct working models of 2.5cm diameter, 
so 2.56cm screen discs were used. 

A screen had to be selected, which had a porosity 
such that it would give the same effect on the 
momentum of the wind as the wind turbine itself. 
The following procedure was used to select the 
screen best suited for the experiment. 

The tv/o working models (28cm diameter) of the 
wind turbines were placed in the Boundary Layer 
Winci Tunnel, with one directly behind the other. 


Both were started, and their rotation speeds were 
“measured, as the distance between them was varied 
from 4 diameters to 32 diameters. The rotational 
speed was measured with a Strobotac. Since the 
Tip Speed Ratio was constant for this particular 
turbine, the rotational speed "w" was directly 
proportional to the wind velocity (Uoo) . 

Seven different types of screens, ranging from a 
very fine mesh to a quite heavy mesh, were then 
substituted for the upstream windmill. 

The distance between the screen and the downstream 
wind turbine was then varied from 4 to 32 diame- 
ters. At each position, the rotational speed of 
the down wind turbine measured. It was then 
possible to select the screen which modelled the 
wind turbine drag, and this screen had a porosity 
of 0.48, Twenty screens, each with a diameter 
of 2.5cm and troposkein shape, were then manufac- 
tured , 

The drag ratio of the screen was measured using 
a small balance that was constructed and calibra- 
ted using a solid disc which had a known drag 
coefficient of 1.17. The drag coefficient for 
the wind turbine was estimated from the power 
coefficient value, and was found to be about .8 
for the range in which the windmill was operating. 
The screen was then mounted on the balance, and 
its coefficient of drag determined to be approxi- 
mately .75. 

The correct screen had now been selected, and it 
had been established by two methods that this 
screen had the same effect on the wind momentum 
as a wind turbine. Therefore, if the 2,5cm 
screen discs were placed in an array, the effect 
on the wind velocity should be similar to the 
effect caused by an array of wind turbines. 

Four sheets of masonite were placed in the tunnel 
and secured to the floor, and the array of screens 
were then placed in the Boundary Layer Wind 
Tunnel. The masonite had holes drilled every 
ten centimeters in each direction, so that the 
size of the array could be varied from 4 diame- 
ters to 36 diameters. In all cases, measurements 
were obtained using an array consisting of five 
rows perpendicular to the wind flow, and three 
parallel to the wind flow. 

Velocity values were recorded downstream from the 
array as the distance between adjacent windmills 
was changed from 4 to 32 diameters. The velocity 
was measured at the location where th_e sixth row 
of the array would have been placed. The vertical 
height of the velocity measurements was equal to 
the mean height of the porous screen models. 
Measurements were taken every ten centimeters 
along the horizontal, by using transversing gear 
located in the Boundary Layer Wind Tunnel, By 
entering the required information into the compu- 
ter, the data was automatically obtained for a 
sample time of 60 seconds, A hot wire anemometer 
was used to obtain the values of the velocity^ 
which were compared to the values of velocity 
when no models were present. The value for no 
models present was obtained from the velocity 
profile which was taken at each location before 
the array was placed in the wind tunnel. The 
fact that the profiles taken at each location 
were similar, was also an indication of the fact 
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that the turbulent boundary layer did not deter- 
iorate over the masonite. 

DISCUSSION OF WAKES 

In the analysis of the wake downstream of a 
Darrieus rotor, two distinct approaches may be 
used. In the near field, the decay may be con- 
sidered to be that of an interacting discreet 
vortex array. For a single two-dimensional vis- 
cous vortex, the circulation (F) decays accord- 
ing to the following equation: 

r = [l-EXP (-r^/4Vt)] 

and this would infer that for the vortex wake 
from a turbine, ' time* was more important in the 
decay process than distance, although both are 
related by the convection velocity of the vor- 
tices, Batchelor (1964) discussed the axial flow 
in trailing line vortices with application to 
trailing vortices from lifting surfaces. An im- 
portant result of this study showed that the 
axial velocity defect varied as : 

{Uo=^")n«x - ^ log[v) 

Where Co is (2 tt) ^ times the vortex circulation 
at large values of radius. For a turbulent wake 
downstream of a circular disc, using the mixing 
length theory (see Schlichting, 1968) , the velo- 
city deficit varies as: 

(U-u) - x ’^/3 

Grainger (1966) studied three-dimensional vortex 
flow and approximated the vorticity by an infi- 
nite power series of the form: 

n(r,x) = f(x) 

An interesting solution to the steady flow vorti- 
city equation was then obtained, and in particu- 
lar, the variation of the vorticity at the centre 
line of the vortex was shown to vary in the form 
of the elliptical function of the Weierstrass 
canonical form. One disadvantage of this type of 
solution when applied to the study of the decay 
of a vortex sheet shed from a Darrieus rotor, is 
that the mutual vortex effects which could also 
cause vorticity amplification, have been ignored. 
Another effect also excluded in these studies, 
which is important when considering the turbine 
wake decay, is the ground effect. At large dis- 
tances downstream of the turbine, the vortex wake 
approaches an established turbulence flow state, 
Sforza (1970), Counihan (1971) and Lemberg (1973), 
considered turbulent wall wakes behind various 
shaped obstacles. The main conclusion of these 
studies, was that for a turbulent wake, the centre 
line velocity deficit varied inversely as the dis- 
tance downstream of the turbine (X"^). These 
established results for the decay of wakes, will 
be used for comparison with measured decay rates 
in the following section. 

RESULTS 
Near Wake 

A vortex model describing the flow field about a 
Darrieus rotor was introduced. In the model, the 
flow field was assumed to consist of a regular 
intermittent coherent pattern of finite vortex 


elements, shed from the rotating blades every 
half a revolution. The circulation strengths of 
the vortices in the vortex model depended on the 
turbine parameters, such as Tip Speed Ratio, so- 
lidity, aerofoil blade section, and blade posi- 
tion relative to the incident airstream. The 
vortex model also satisfied Kelvin's theorem con- 
cerning continuity of vortices. Figures 2 and 3 
show typical vortex arrays used in the model. 

In this initial study, it was shown that for a 
two-bladed Darrieus rotor, the vortices shed from 
the 'equatorial* element of the upstream blade 
intersected the downstream blade after 150 of 
rotation, when the rotor Tip Speed Ratio had a 
value of IT, Further computer studies on this 
vortex model of the Darrieus rotor are planned 
in the future. 

Far Wake 

The first tests were to determine the correct 
screen with a given porosity, to model the 
Darrieus rotor operating at the design Tip Speed 
Ratio of 5,5, Two similar turbines were placed 
in a line parallel to the ambient air flow, so 
that the downstream turbine was affected by the 
wake of the upstream turbine. The downstream 
turbine was then moved downstream in increments 
of four dicimeters, and the various rates of ro- 
tation of the downstream turbine recorded. The 
forward turbine was then replaced by troposkein 
shaped screens of different porosity, and the 
rate of rotation was again recorded for the down- 
stream wind turbine. Figure 5 shows the test 
results for the rate of rotation of the down- 
stream turbine when initially a Darrieus turbine 
was upstream, creating blockage, and also with a 
screen with a porosity of 0.48, which was used 
in subsequent experiments placed upstream. 

Tests were now made on an array of ’wind tur- 
bines* modelled by one inch diameter porous 
screens. Figure 6 shows the (3 x 5) array of 
fifteen screens in a square pattern array. The 
velocity measurements were made, as shown, in 
the sixth row. Figure 7 shows a typical wake 
traverse along the sixth row, when the screen 
array spacing was eight diameters. The maximum 
velocity deficit was taken as a measure of the 
wake effect. For example, when the square 
spacing was eight diameters and the free stream 
velocity 18.9 ft/sec., the velocity deficit was 
4.4 ft/sec. This implied a percentage velocity 
deficit based on the free stream velocity of 
23.3%, and a power deficit of 54.8%, Figures 8 
and 9 show the estimated velocity and power 
deficits respectively, due to the wake effects 
of the five rows of screens, representing an 
array of Darrieus turbines for various square 
spacing, with the measurements being made in the 
sixth row position. Also in Figure 8, the wake 
decay is compared with the wake decay predicted 
by the various theories discussed in the pre- 
vious section. 

CONCLUSIONS 

The theory and experimental measurements on the 
wake of a Darrieus turbine and an array of tur- 
bines were studied: In general, the development 

of the turbine wake may be classified into 3 
regions: The first region consisted of the for- 
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mation and rolling up of the concentrated vortex 
__ , sheet f and ^en the decay by molecular tur- 

buisnt diffusion. Finally, a homogeneous, fully 

developed turbulent state was achieved far down- 
stream of the turbine. For a two-bladed Darrieus 
rotor or straight-bladed vertical axis turbine, 
when the vertical plane of the rotor was normal 
to the flow direction, the lift force on the 
blades was zero. At any other angle, the local 
irictdence of the aerofoil blades at the ^equator* 
for instance, could be calculated by the follow- 
ing expression! 

Z (l-a TS R I Sine I ) Sin9 

- (1-0 Jiin6|)Cos6 + Tsr 

-- --- Ftir the time per io d for the rotor to rotate half 
- “a revolution, a closed vortex array was shown to 
be shed from the aerofoil blades, and was convec- 
ted downstream, changing shape due to the mutual 
vortex interaction and the effect of the ambient 
air flew. 

In the experimental study of the wake far down- 
stream of an array of Darrieus wind turbines, 
represented by troposkein shaped porous screens, 
the results indicated that the ’turbines* should 
be spaced at least forty diameters apart to 
avoid mutual power depreciation greater than ten 
per cent. This predicted value of suitable 
spaaing of turbines, compares favourably with 
the analytical estimates of Templin (1974) , 
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NOMENCLATURE 

a interference or blockage factor for vertical 
axis wind turbine 

B number of windmill blades 
c blade chord 

Co vortex circulation times (27T) ^ 

Cl blade section lift coefficient 
d diameter of wind turbine 

r radius of a point in a vortex from the centre 
line or radius of a point on a turbine 

h height of blade element of wind turbine 
t time 

TgR Tip Speed Ratio (HR/Uco) 

ut total local airspeed relative to blade 

Uoo ambient (approaching) wind velocity 

a angle of attack 

X distance downstream of turbine 

6 angle of blade rotation for vertical shaft 
windmill 

V kinematic viscosity 

p density of fluid 

0 solidity (BC/2R) 

Q blade rotational speed or vorticity 
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r circulation 
To initial circulation 
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Nn RECTANGULAR VORTEX SYST 
A SINGLE ROTATING AEROFOIL 
IN A CROSS FLOW 








RGURE 4 SCHEMATIC VIEW OF T\WO 

BLADED DARRIEUS ROTOR 
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FIGURE 6 


Measurements made along this line 



SCHEtlATIC OF 'WINDTURBINE' SCREEN ARRAY 
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MEASUREMENTS MADE 'IN SIXTH ROW 





MEASUREMENTS MADE IN SIXTH DOWNSTREAM 


QUESTIONS AND ANSWERS 
T.E. Base 


From: T,A. Egolf 

Q: Isn’t it true that if you reduce your time step interval to a sufficiently small 

At, that the unsteadiness due to the blade impact with the shed vortex sheet is 
significantly reduced, and if unsteady aerodynamics were used, the effect would 
be further reduced? 

A: If the time increment is too shorty the shed vortices from the blade will not he 

formed into a coherent structure , 


From: T, Sullivan 

Q: How do you determine wake position for the near wake model? Are any measurements 

of this wake position available? 

A' The wake ie shed from the wing and the subsequent position is determined from the 
induced flow field of the other vortex elements and are free stream, i.e. 




No* Vortices 
S 

i = l 


M . -h U 


I do not know of any published detailed near wake velocity traces or vorticity dis 
trihutions or flow wsual'Lzation studies » 
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